On positivity of principal minors of bivariate Bezier 

collocation matrix 



Gasper Jaklic a , Tadej Kanduc*' b 

a FMF and IMFM, University of Ljubljana and PINT, University of Primorska, 
Jadranska 19, Ljubljana, Slovenia 
b Turboinstitut d.d., Rovsnikova 7, Ljubljana, Slovenia 



Abstract 

It is well known that the bivariate polynomial interpolation problem at do- 
main points of a triangle is correct. Thus the corresponding interpolation 
matrix M is nonsingular. L.L. Schumaker stated the conjecture, that the 
determinant of M is positive. Furthermore, all its principal minors are con- 
jectured to be positive, too. This result would solve the constrained inter- 
polation problem. In this paper, the basic conjecture for the matrix M, the 
conjecture on minors of polynomials for degree < 17 and for some particular 
configurations of domain points are confirmed. 
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1. Introduction 

Positivity of determinants (or minors) of collocation matrices is an im- 
portant property in approximation theory. Nonsingularity of such a matrix 
implies existence and uniqueness of the solution of the associated interpola- 
tion problem. Positivity of principal minors, or even total positivity is used 
in the proofs of some well-known results, see [l], 0, e.g. 

Recently, nonsingularity and principal minors of such matrices for poly- 
nomial interpolation at Padua-like points and for interpolation by bivariate 
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Bezier patches were studied in |4| and |5|]. A related problem is a construction 
or an approximation of Fekete points for a given domain, i.e., the interpo- 
lation points, which yield the maximal absolute value of the Vandermonde 
determinant ((J. 

It is a well known fact that bivariate Bernstein polynomials {Bf- k \i + jj r k=d 
of degree d form a basis of the space of bivariate polynomials of degree < d. 
Let l d — {{i,j,k) : i+j + k = d, i, j, k e N U {0}}. For every ( d + 2 ) points 
in a domain, which do not lie on an algebraic hyper-surface of degree < d, 
the corresponding interpolation problem is correct (0]). In particular, the 
problem is correct for domain points 

T>d,T ■= {(i/d,j/d,k/d) : (i,j,k) e l d }, 

expressed in barycentric coordinates with respect to a given triangle T. Thus 
the corresponding interpolation matrix M := [-B^(£)ke£>d T,v^d * s nonsingu- 
lar. In [8j, a theorem was stated that by choosing an arbitrary nonempty 
subset J C T>d,T and the corresponding set of indices T, the sub matrix Mr := 
[B%(£)]££j t p£r is nonsingular for all d < 7, and furthermore, det Mr > 0. The 
authors of [8] verified the theorem by computer and provided a proof of non- 
singularity of principal matrices for some special configurations of domain 
points. For more details see [5(. 

A straightforward way of verifying the conjecture by computing principal 
minors of M is time consuming due to the exponential growth of the number 
of subsets that need to be analysed, and cannot be done for d > 7 using 
current computational facilities. 

Matrices with positive principal minors are known as P-matrices. Much 
on them is known, see 0, H3], e.g. Unfortunately, those tools could not be 
applied for the study of the problem at hand, so a different approach will be 
used. 

In this paper, the conjecture on positivity of determinant of the bivariate 
Bezier collocation matrix M is confirmed. Furthermore, the result is proven 
for an arbitrary T C for d < 17. Thus the constrained Lagrange inter- 
polation problem has an unique solution. This covers all the cases useful in 
practice, since it is well known that Lagrange polynomial interpolants of high 
degrees have undesired properties. Some particular configurations of domain 
points are analysed. A conjecture for an exact lower bound of det Mr is 
stated. 
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2. Main results 



Let i be a weak 3-composition of an integer d, i.e., i = k), such that 
|i| := i+j + k = d and k G NU{0}. Let X d := {i}^ =d be a set of all weak 
3-compositions of the integer d. The set X d consists of ( d ~^ 2 ) compositions. 

Let T be a triangle in the plane P. Every point v G P can be written in 
barycentric coordinates v = (u, v , w) , u + v + w = 1, with respect to T. The 
Bernstein basis polynomials of total degree d in barycentric coordinates are 
defined as 



Bf(v) :=Bt lk (u,v,w):- 



d\ 



v 



i\j\k\ 



-u l v j w k , 



d. 



Here the standard multi-index notation and a convention 0° = 1 are used. 

(z) 

Let us denote the subset of all compositions with z zeros by X d C I d , 
z = 0, 1, 2. Let £i := := i/d be a domain point of the triangle T, repre- 
sented in barycentric coordinates. A domain point & is boundary if at least 
one of its barycentric coordinates is zero, i.e., i G T d UX d . 

We are now ready to present the conjecture, stated in that will be 
tackled in this paper. 

Conjecture 1 ([5]). For a given triangle T and every nonempty set T = 
{ii, i^i ■ ■ ■ , i n } C X d , the matrix 



Mr 



[B}{ti)\ 



Bi 2 (£i 2 ) 



BtHh 

Ba (62 



Bi (Cin , 



is nonsingular. Furthermore, detMp > 0. 

A confirmation of Conjecture [1] would imply the following. Let V d X d 
and let £({Bf}i e r) be the given interpolation space. Then the interpolation 
problem for the points {£i}ier in the domain would be correct, i.e., the 
interpolant would exist and be unique. An example is shown in Fig. [H The 
conjecture is important for interpolation with spline functions, since some 
degrees of freedom are determined through the smoothness conditions and 
the rest through the interpolation conditions (see 11, [Ej, e.g.). 
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Figure 1: Let boundary control points Cj) : iel] Ulj , q eK}ofa Bezier surface 
p and real values z 2 n, Zm, ^112 be fixed. Conjecture Q] ensures that interior control points 
{(£i> c i) : i S , Ci 6 M} can be set in such a way that p(£i) = i G X4 . 



Note that the matrix Mr is not symmetric. The determinant of M-p is 
independent of the ordering of elements of F as long as the same ordering for 
rows and columns is used. It is common to use the counter-lexicographical 
ordering K-iex, 

(tZ, 0,0), (d- 1,1,0), (d- 1,0,1), (d-2,2,0),..., (0,0, d), 

but a particular ordering of elements in X^, which yields a block lower trian- 
gular matrix Mx d (see 12] ), will be more convenient. A linear ordering is 
defined as: i >~b j if one of the following holds true: 

1. i G l ( d Zl) and j G X ( d Z2) for z x > z 2 , z t G {0, 1, 2}, 

2. i, j G X d z) for z G {0, 1, 2} and sgn(t) K-iex sgn(j'), 

3. i, j G xj° for z G {0, 1, 2} and sgn(i) = sgn(j) and i K-icx j. 

Here sgn(i) = sgn(i,j, k) := (sgn(i), sgn(j), sgn(fc)). The ordering implies 
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that the matrix Mj d has a structure 



M t(2) 



* 



Mux) 



d 



* 



* M (o) 



■d 



where 



M {{d,0,0)} 



Mum 




d 




Mux } 



d 



and e := {£ ii0 G X d : i, j > 1}. 

Thus the problem of verifying the positivity of principal minors of the 
matrix Mj d is reduced to each diagonal block matrix separately. The matrix 
M e is an univariate Bezier collocation matrix and by [13| it is totally non- 
negative with positive principal minors. Therefore, the problem reduces to 
the study of a ( d ~ 2 1 ) x ( d 2 l ) matrix Muo), which, unfortunately, represents a 

d 

very large part of the matrix Mj d for a large d. 

Let us simplify the considered matrix. Let us construct a matrix from 
Mr in the following way: 

• for every column, divide each element, that corresponds to a polynomial 



Clearly, the matrix iVr is a principal submatrix of Nx d - Since sgn det Mr = 
sgndetiVr, Conjecture [U holds for Mr iff it holds for iVr- 




det N r 




det M r . 



(1) 
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The matrix has some nice properties. It consists only of non-negative 
integers, thus determinant computations are exact. The matrix Nx d has 
a simpler structure than Mr d and is closely related to combinatorial ob- 
jects. Therefore some properties of the matrix Mr d will be proven via Nx d . 
Note that some of the properties are not preserved by the transformation 
M r JV r . 

As an example, diagonal blocks of the matrix Mj 4 and the corresponding 
matrix iVj 4 are shown in Table [TJ 



rc/ 4 



Mr 



Nr 



{(4,0,0)} 
{(0,4,0)} 
{(0,0,4)} 
{(3,1,0),(2,2,0),(1, 
{(3, 0,1), (2,0, 2), (1, 
{(0, 3,1), (0,2, 2), (0, 
{(2,1,1)} 
{(1,2,1)} 
{(1,1,2)} 



3,0)} 
0,3)} 
1,3)} 



i 

128 



32 



1 

54 27 6 
32 48 32 
27 54 
3 3 
6 3 
3 6 



6 



6 
3 
3 



256 

27 9 3 

16 16 16 

3 9 27 

4 2 2 

2 4 2 

2 2 4 



Table 1: Diagonal blocks of the matrices Mx 4 and Nx 4 by using the ordering >-(,. 



Now we are ready to present one of the main results of the paper. 

Theorem 2. Let d < 17. Then Conjecture^ holds true, i.e., detMr > for 
every nonempty subset V C Id- 
Proof. Fix d, 1 < d < 16. The matrix Mx d + Mj is symmetric and positive 
definite since its Cholesky decomposition exists. Therefore the matrix Mr d 
is positive definite, i.e., x T Mx d x > for all x ^ 0. Thus all principal 
submatrices of M% d are positive definite too and all principal minors are 
positive. 

For d — 17, the matrix Mx d + M-£ has three negative eigenvalues. The 
matrices M T m and M (i) are P-matrices, thus the problem reduces to the 
study of the matrix M (o). Since the Cholesky decomposition of the matrix 
M (o) + M T (0) exists, the rest of the proof is similar to the first part. □ 

d I'd 
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For d = 18, the smallest eigenvalue of Mj d + M£ is approximately —1.1 x 
10~ 7 . For d > 18, the number of negative eigenvalues of the matrix Mj d +Mj 
increases with d. Therefore, this approach can not be used to prove the 
conjecture in general. However, our result covers all the cases, important 
in practice, since only Lagrange polynomial interpolants of low degrees are 
useful. 

Theorem |5] implies the following result. 

Theorem 3. Let T C 1^ and let d < 17. Then for any {zi}i e r, there is an 
unique polynomial of the form 

ier 

such that 

p(&) = Zi, ie r. 

Remark 1. Theorem[3]generalizes 0, Thm. 3]. Its proof avoids computation 
of all subdeterminants, as was the case in [§]. 

Remark 2. Note that the matrix Nj d is positive definite only for d < 4. 

Remark 3. The approach for proving Theorem [2] can also be used to prove a 
similar result for the trivariate case: the trivariate Bezier collocation matrix 
for d < 16 is positive definite, and thus the constrained polynomial Lagrange 
interpolation problem is unisolvent. 

Now let us consider some particular configurations of domain points (and 
the corresponding choices of T) for d arbitrary. 

Theorem 4. Let d be arbitrary and let F satisfy one of the following as- 
sumptions: 

1. |r| < 2, 

2. let one of the components of k) be fixed for all elements in T, 

3. r = i d , 

4. Y = Tf\ 

5. rcxfuxf, 

6. r = Ti ur 2 , where Ti is one of the sets, defined in 1., 2. or 4-, and T 2 
is a set in 5. 
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Then detM r > 0. 



Proof. Let |T| < 2. For |T| = 1, the matrix Mr is a positive number. 

Now let T = i 2 }. Since the largest element of every column in M-p is 
on the diagonal of Mr, 



detMr 



> o. 



Let one of the components of (i, j, k) be fixed. Without loss of generality 
we may assume that it> = k), %g + ji + k = d, £ G {1,2,..., |T|}. By 

dividing each element of N-p by k k and multiplying each column by a proper 
constant, the matrix Np transforms to an univariate Bezier collocation ma- 



trix, which is a P-matrix by [13 



Now let us consider the case T = X^. It follows from 12J and (|T]) that 



(,\ min{d,3} / d—k+1 \ (fc) 

•) n n ^ (n ^ +i)( "- i) >o. 

' k=l \ 1=1 J 



(2) 



Positivity of det M T (o> follows from 



det M Xw = det M t(2) det M T m det M t(0) > 0. 

x d x d x d 

For T C 2^ U2j and r = Ti U T 2 , the result follows straightforwardly 
from the structure of the matrix Mr. □ 

Remark 4. The set T in Theorem H(2) corresponds to domain points in the 
triangle T, lying on a line parallel to some edge of T. 

The nonsingularity of the matrix Mj d follows from the fact that Bernstein 
polynomials form a basis of the space of bivariate polynomials of degree < d. 



Although the proof that det Mj d > may seem easy, most of the paper [12 
is dedicated to the derivation of determinant formula in a closed form (j2]). 

The subset of decompositions T in Theorem H^5) corresponds to the in 
terpolation problem at boundary domain points of the triangle T. 
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3. Concluding remarks 

In the rest of the paper we will state some conjectures which expand the 
Conjectured! and examine some interesting properties of considered matrices. 



Conjecture 5. For d fixed, 



mm detM r = det M Trl , 

rci d 

ry0 



where det Mx d is given in ([2]) . 
Conjecture 6. For £ e N let 



Then 



£ u , d = 3£ 
n d := { (£+l) i+1 £ 2e , d = 3£+l . 
(£+l) 2e+2 £ e , d = 3£ + 2 



min det Nr = rid- 

rcx d 
r^0 



Conjectures O and were verified by a computer for d < 7. Let us prove 
the latter conjecture for |T| < 2 and arbitrary d. We will need the following 
lemma. 

Lemma 7. Let x = (x, y, z) e R 3 and fix i = k) G X<j. Let the function 
f{x) := x l be defined on 

Q = {(x,y,z) : x + y + z = d, < x < d, < y < d — x}. (3) 

Then f has an unique maximum at i and 

max fix) = i 1 . 

Proof. Let T = ((0,0), (d, 0), (0, d)) be a triangle in the domain and let us 
define Bernstein polynomial Bf on T. Since Bf = (i)/d d f and Bf has an 
unique maximum in the barycentric point £i, the proof is complete. □ 

Proposition 8. Conjecture^ holds true for \Y\ < 2. 
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Proof. Let |T| = 1 and let g(x) := g(x,y,z) := x x be a function, defined on 
the domain Q as in ([3]). We are looking for 

a := min g(x). 

xennz 3 

An unique local minimum of g in the interior of Q is obtained as a so- 
lution of the normal system dg/dx = 0, dg/dy = 0, and it is reached at 
(d/3,d/3,d/3). This is a global minimum since g(d/3,d/3,d/3) < g(x) for 
all x at the boundary of Q. 

If d = (mod 3), then a = {d/3) d . 

Let us examine the case d = l(mod 3). Then + for 

£ G N U {0}. By the symmetry of the function g and since at least one 
component of x G Q fl Z 3 is greater or equal to £ + 1, it is enough to consider 
the case £ + 1 < x only 

Let us define 

Q x := {x = (x, y, z) G Vt : £ + 1 < x}. 

Since g has no extreme point in Q x , the minimum value is reached at the 
boundary of Q x . Then the minimum is a = (£ + l) e+1 £ 2i and it is achieved 
at + 1, £, £). For £ + 1 < y and £ + 1 < z, the derivation is analogous. 

The case d = 2 (mod 3) is similar to the previous one. Since N^y = g(i), 
i G Td, and rid = a, the conjecture for |T| = 1 is proven. 



Now let us consider the case |T| = 2. Let us show that 

it 1 = detN {il} < detN {iui2} = i^€? - i* 1 ** 2 

for every i x , i 2 G l d , i\ ^ *2- 

By Lemma [7] it follows that i^ 1 < i* 1 and ** 2 < i % 2 2 — 1, thus 

< t{ 1 (*« - 1) - i* 1 ** 3 = detiV {ilii2} - detiV {il} . 

□ 

Let us denote an eigenvalue of the matrix Mx d with index i and multi- 
plicity n by \ff^ ■ An interesting interlacing of eigenvalues can be observed. 

Conjecture 9. The spectrum of the matrix Mj d consists of d different real 
eigenvalues 

1 \<d> \<d> . \<ci> \<d> ^_ ^_ \<d> n 

i - A l,3 > A 2,3 > A 3,4 > A 4,5 > ' ' ' > A d, d+1 > U - 
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0.0 0.2 0.4 0.6 0.8 1.0 

Eigenvalue 



1.2 



Figure 2: Eigenvalue interlacing for M% d , d < 20. 



The eigenvalues satisfy a criterion, similar as in Cauchy interlacing theorem, 

A<* d+1> > A<f > Ag£> i = l,2,...,d. 

The conjecture was verified by a computer for d < 20. Note that positivity 
of eigenvalues would follow, if one would prove that Mx d is a P-matrix. 

As an example, the spectra of Mj 10 and Mj 10 + Mj io are shown in Table [2] 
and El respectively. The interlacing property for d < 20 can be observed in 
Fig. El 
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Eigenvalue 


Multiplicity 


1 


3 


0.9000 


3 


0.7200 


4 


0.5040 


5 


0.3024 


6 


0.1512 


7 


0.06048 


8 


0.01814 


9 


3.629 x 10~ 3 


10 


3.629 x 10" 4 
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Table 2: Spectrum of the matrix Mx 10 . 



Eigenvalue 


Mult. 


Eigenvalue 


Mult. 


Eigenvalue 


Mult. 


2.341 


1 


0.2854 


1 


7.224 xlO" 3 


1 


2.308 


2 


0.2675 


1 


7.038 xl0~ 3 


2 


1.856 


2 


0.2585 


1 


6.565 xl0~ 3 


1 


1.841 


1 


0.2568 


2 


6.560 xl0~ 3 


1 


1.475 


1 


0.1225 


1 


6.399 xl0~ 3 


2 


1.412 


2 


0.1179 


2 


6.220 xl0~ 3 


2 


1.303 


1 


0.1076 


2 


6.209 xl0~ 3 


1 


1.001 


2 


0.1033 


2 


6.752 xl0~ 4 


2 


0.9377 


1 


0.1028 


1 


6.600 xlO^ 4 


2 


0.8855 


2 


0.03647 


2 


6.574 xlO~ 4 


1 


0.6082 


1 


0.03531 


1 


6.160 xlO~ 4 


2 


0.5742 


2 


0.03253 


2 


6.079 xlO~ 4 


1 


0.5221 


2 


0.03241 


1 


5.951 xlO~ 4 


1 


0.5180 


1 


0.03126 


1 


5.923 xlO~ 4 


1 


0.2995 


2 


0.03115 


2 







Table 3: Spectrum of the matrix M% w + Mj io together with eigenvalue multiplicities. 



12 



References 

[1] C. de Boor, A practical guide to splines, Springer, 1994. 

[2] M. Gasca, Spline functions and total positivity, Revista Mathematica 9, 
1996. 

[3] M. Gasca, C.A. Micchelli, Total positivity and its applications, Springer, 
1996. 

[4] L. Bos, S. De Marchi, S. Waldron, On the Vandermonde determinant of 
Padua-like points, Dolomites Research Notes on Approximation, vol. 2 
(2009), 1-15. 

[5] M.-J. Lai, L.L. Schumaker, Spline functions on triangulations, Cam- 
bridge University Press, 2007. 

[6] L. Bos, S. De Marchi, A. Sommariva, M. Vianello, Computing multivari- 
ate Fekete and Leja points by numerical linear algebra, SIAM J. Numer. 
Anal. 48 (2010), 1984-1999. 

[7] G. Jaklic, J. Kozak, M. Krajnc, V. Vitrih, E. Zagar. Three-pencil lattices 
on triangulations, Numer. Algorithms 45 (2007), 49-60. 

[8] G. Nurnberger, V. Rayevskaya, L.L. Schumaker, F. Zeilfelder, Local La- 
grange interpolation with C 2 splines of degree seven on triangulations, in 
Advances in Constructive Approximation: Vanderbilt 2003, M. Neamtu 
and E. Saff (eds.), Nashboro Press, Brentwood, (2004) 345-370. 

[9] K.G. Murty, On a characterization of P-matrices, SIAM J. Applied 
Mathematics, 20 (3), (1971), 378-384. 

[10] G.E. Coxson, The P-matrix problem is co-NP-complete, Math. Pro- 
gramming, 64(2) (1994), 173-178. 

[11] G. Jaklic, J. Modic, On properties of cell matrices, Appl. Math. Corn- 
put., 216 (2010), 2016-2023. 

[12] J.M. Brunat, A. Montes, The power-compositions determinant and its 
application to global optimization, SIAM J. Matrix Anal. Appl., 23(2) 
(2001), 459-471. 



13 



C. de Boor, R. deVore, A geometric proof of total positivity for spline 
interpolation, Math. Comp. 45(172) (1985), 497-504. 



14 



